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Abstract. We propose that the mushroom-shaped structure of the Galactic worm GW 123.4-1.5 is created by a 
cloud collision with the Galactic gas disk. A hydrodynamic simulation shows that a mushroom-shaped structure is 
created after the cloud crosses the Galactic midplane. The lifetime of the mushroom-shaped structure is of order 
the dynamical time scale of the disk, ~ 10 7 years. We find that the velocities across the cap of the mushroom- 
shaped structure in the simulation are consistent with the observed values. The simulation also predicts a structure 
on the opposite side of the Galactic plane which is created by the Kelvin-Helmholtz instability after the cloud 
passes through the disk. 
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1. Introduction 

The Canadian Galactic Plane Survey (Taylor et al. 2003) 
revealed that the Galactic worm candidate GW 123.4-1.5 
was an unusual mushroom-shaped cloud (English et al. 
2000). It is hundreds of parsecs in size and unrelated to 
conventional shell or chimney structure. The mass of the 
cap is about 4 times greater than that of the stem, and 
the total mass is estimated to be about 1.55 x 10 5 Mq, al- 
though such an estimate depends strongly on the distance 
in which there is some ambiguity. A position- velocity map 
in the upper portion of the cap shows no velocity gradi- 
ent. However, the lower portion of the two lobes of the cap 
that extend back toward the Galactic plane are blueshifted 
with respect to the central cap region by 5 km/s. 

One possible origin for the mushroom-shaped cloud is 
the rise of buoyant gas. This model was studied by English 
et al. (2000), in the form of a supernova event, and Avillez 
& Mac Low (2001), in the form of bubbles rising from 
hot gas reservoirs. They showed that a mushroom-shaped 
cloud was created by the sweeping effect of the buoyant 
bubble when the bubble rises from the Galactic plane. The 
gas swept up by the bubble cooled down and increased in 
density. 

In this paper, we propose another possible origin of 
the mushroom-shaped structure: a cloud collision with the 
Galactic disk. This scenario was first studied by Tenorio- 
Tagle et al. (1986, 1987). However, no one has yet per- 
formed a simulation demonstrating that a cloud impact 
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would generate the mushroom-shaped structure like GW 
123.4-15. 



2. The modeling of a gas collision with the 
galactic disk 

We perform a hydrodynamic numerical simulation. As a 
first step, we simply assume adiabatic gas and axial sym- 
metry. We model a local region of the Galaxy and use the 
axially symmetric cylindrical coordinate (R,z), where R 
represents the distance from the axis and z the vertical 
distance from the midplane of the Galaxy. 

According to observations, the vertical gravitational 
field of the Galaxy is g z = —az near the midplane, where 
a ~ 2/3 x 10 -29 s -2 for the solar neighborhood (Spitzer 
1978; Bahcall 1984; Tenorio-Tagle et al. 1987). The hydro- 
static equilibrium density of the interstellar isothermal gas 
near the midplane is, 



p(z) = p cxp [-(z/H y 



(1) 



where po is the density at the midplane. In equation l|T)l. 
the scale length is 



H 



(2) 



where 7 is the specific heat ratio, assumed to be 5/3, and 
c s o is the sound speed at the midplane. Since the gravita- 
tional field of the Galaxy gradually tends to be constant 
when z is large, we simply let 



-aHg tanh(z/iJ g ), 



(3) 
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which tends to g z 



-az when z <C H g and to g z 



-aHg 

when z ^> H g . We take H g — l.bHo in the simulation. 
Some calculations were done by using larger values of H g , 
but the results were not sensitive to it. 

We also assume that the interstellar gas has two tem- 
peratures at the initial time so that 



T(z) = T + 0.5(71 - T ) 



1 + tanh 



zt 



(4) 



where T is the temperature as a function of z. The temper- 
ature at the midplanc is the constant To. It becomes an- 
other constant T\ at \z\ = Zt with a transition length of Zd- 
We take T x = 10T , z t = 0.5H„ and z d = Q.IH Q . Although 
observations currently indicate a multi-temperature galac- 
tic atmosphere (e.g., Lockman & Gehman 1991), a two 
temperature model of the Galactic disk is a first step that 
adds a level of sophistication beyond the isothermal at- 
mospheres explored by Tenorio-Tagle et al. (1986; 1987). 

As an initial condition of the Galactic plane, we as- 
sume a hydrostatic equilibrium of the gas determined by 
the gravity of equation and the temperature of equa- 
tion |0J. When z <C H g , the gravitational field is propor- 
tional to z and the density has a Gaussian distribution 
like equation There is a transition region of the den- 
sity at zt = 0.5-f/o due to the two temperature model 
of equation |@J. Because the gravity tends to a constant 
where z 3> H g , the density structure finally tends to an 
exponential atmosphere at large z. The hydrostatic equi- 
librium of the two temperature atmosphere makes a cold 
dense sheet which is sandwiched between hot gas layers. 
As we will see later, this dense sheet plays an important 
role in making a stem structure of the mushroom. 

The density distribution in our simulation is not ex- 
actly the same as that shown by Dickey & Lockman (1990) 
who found that the z-distribution of the disk density was 
approximated by a summation of two Gaussians which 
have different velocity dispersions. The distribution they 
found does not show a clear discontinuity in the temper- 
ature. However, their distribution is the average density 
along the line of sight. In this paper, we assume that the 
cold gas matter is compact and can be locally well sepa- 
rated from the warm gas. In addition to this, we assume 
that the cold gas is located in a sheet near the midplane 
of the disk. If the cold gas clouds along the line of sight 
have a distribution of heights about the midplane of the 
Galaxy, then a transition region of density or temperature 
would become unclear in the observations. 

We also note that the temperatures (To an d Ti) of 
the gas in this simulation are effective temperatures. 
Observations show that the scale height of the Galactic 
interstellar medium is much larger than the scale height 
that is estimated from the temperature of the gas. The 
effective temperature that makes the large scale height is 
considered to be caused by turbulent motions. 

Into this equilibrium atmosphere, we introduce a cloud 
with density (p c ) and z-component of velocity (v zc ): 



p c (R,z) = PcQ 



1 



— tanh 
2 



r c (R,z) - r u 
rd 



v zc {R, z) = v zc0 



where 



2 \ r d 



r c (R,z) = ^R 2 + {z- zc) 2 



(6) 



(7) 



(5) 



and z c is the central position of the cloud on the axis. In 
equations 10 and ©, p c o and v zc o are constants which pa- 
rameterize the initial cloud density and the z-componcnt 
of the cloud initial velocity. The parameters r wc and r wv 
determine radii within the cloud for the density and veloc- 
ity distributions respectively, and is a transition length. 
We use z c = —2.5H , r wc = 0AH , r d — O.IH and 
r wv — r wc + 2r d . We let r wv be a little larger than r wc in 
order to avoid a velocity gradient within the dense cloud. 
(If we take r wv — r wc , a part of the cloud near the edge 
moves slower than the center of the cloud.) 

In this simulation, we take Hq as a unit of length, c s o 
as a unit of velocity and po as a unit of density. If we use 
the value of a for the solar neighborhood and c s q — 10 
km/s (which corresponds to To ~ 10 4 K) as the turbulent 
velocity in the Galaxy, then Ho ~ 140 pc. The unit of 
time is to = Hq/c s o — 1.4 x 10 year. The important 
parameters in this simulation are the initial density of the 
cloud, p c o, and the initial velocity of the cloud, v zc q. In this 
paper, we will show the result of p c o = po, which means 
that the initial density of the cloud is the same as that 
of the cold dense sheet on the midplane of the Galaxy 
(e.g., ~ 1 cm -3 ). We will show two cases of the initial 
velocity; v zc o — 10c s o ~ 100 km/s and v zc o — 5c s o ~ 50 
km/s. The former is the order of the line-of-site motion 
of a high- velocity cloud (Wakker & van Woerden 1997), 
and the latter corresponds to a intermediate- velocity cloud 
(Danly 1989; Kuntz & Danly 1996). 

We use the CIP method (e.g., Yabe, Xiao & Utsumi 
2001) to solve the finite difference equations of hydrody- 
namics. The minimum grid size is 0.01i/o in both the R 
and z directions. In the i?-direction, we gradually increase 
the grid size where R > 3Hq- The computational region 
is < R < 10H , -5H < z < 10H . The outer bound- 
ary conditions are free boundaries for both the R and z 
directions. 

3. Results 

3.1. The case ofv zc o = 10c sU — 100 km/s 

Figure 1 shows the time evolution of the density. The top- 
left panel shows the initial condition (t = 0). The ini- 
tial atmosphere is hydrostatic equilibrium determined by 
equations © and A gas cloud is input in the atmo- 
sphere according to the equations J^J and ©. 

As the cloud impacts the midplane of the Galaxy, 
it goes through the dense sheet in the Galactic plane 
(t = Q.Qto). After that, the cloud is decelerated by gravity 
(t = 1.45£o) and is broken by eddies of the hot gas flow 
generated behind the cloud after it passes through the disk 
(t > 2.75io). The cloud broken by eddies forms a cap of a 
mushroom-shaped structure. 
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Fig. 1. Time evolution of density in the case of v zc o = 10c s o — 100 km/s . The gray scale shows the logarithmic scale 
of the density which is normalized by the initial density at = 0, po- Arrows represent the velocity vectors which are 
normalized by the initial sound speed at z = 0, c s0 . Velocity vectors are shown only in the region R > 0. 



When the gas cloud passes though the Galactic plane, 
a part of the dense sheet in the Galactic plane goes 
along with the cloud and forms a shell-like structure 
(t = 0.6t ). This shell is elongated and forms the stem 
of the mushroom-shaped structure (t = 1.45t — 2.75io)- 
Moreover, when the hole that was made by the cloud pas- 
sage is refilled with the gas on the plane, a flow appears 
near the axis (t = 1.45io — 2.75to). This flow lifts a part 
of the dense gas in the Galactic plane, which also con- 
tributes to the stem (t = 3.5io)- These processes show the 
importance of the dense sheet in the Galactic atmosphere 
that results from our assumed transition in effective tem- 
perature. During this process, the stem is getting narrow 
(t = 2.75to — 4.5to)- The simulation also shows the de- 
velopment of small structures on the opposite side of the 
Galactic plane (z < 0). The two-sided eddies, evident at 
R = ±H and z = -1.5H in Figure 1 at t = 2.75t , 



are created by the Kelvin-Helmholtz instability after the 
cloud passes through the disk (t = 1.45to)- Finally, the 
cap is diffused by eddies and becomes unclear (t > 4.5to)- 
The lifetime of the mushroom-shaped structure is the or- 
der of the dynamical time scale, t ~ to ~ 10 7 year. The 
movie for the density clearly shows the time evolution. 

Figure 2 shows a snapshot of the column density at 
t = 2.75to- Since this is an axially symmetric simula- 
tion, we can integrate the density perpendicular to the 
R — z plane to get a column density. Here, we have in- 
tegrated the low temperature gas (T < 0.8Ti) when we 
calculate the column density, in order to avoid includ- 
ing the gas which is surrounding the mushroom-shaped 
structure. Sampling the low temperature gas makes the 
mushroom-shaped structure clear to understand, though 
this is not exactly the same column density which is ob- 
tained by observations. The column density of low tern- 
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-3-2-10123 
R/H 

Fig. 2. Column density of low temperature gas [T < 
O.8T1) at t = 2.75t in the case of v zc0 = 10c s0 100 
km/s. The gray scale shows the linear scale of the column 
density normalized by poH . 

perature gas has a clear stem and cap structure which re- 
sembles the mushroom-shaped structure of GW 123.4-15. 
(The movie for the column density also shows its time evo- 
lution.) The width ratio of the cap and stem at this snap- 
shot is about 3:1, and the mass ratio is about 2.4:1. These 
results are consistent with a observational ratio of 3:1, 
and the mass ratio of 4:1. However, the maximum height 
of the mushroom at t = 2.75t is about 5.5H — 770 pc, 
which is greater than the mushroom-shaped structure of 
GW 123.4-15 though the scale depends on the ambiguous 
distance. 

The column density in Figure 2 also shows a smaller 
structure on the opposite side of the Galactic plane 
(around z = —l.5H ). This structure is created by 
the Kelvin-Hclmholtz instability after the cloud passes 
through the disk. The lifetime of this structure is also 
~ 10 7 year. 

Another important comparison with observations is 
the velocity structure of the mushroom-shaped struc- 
ture. The observational velocity, which is obtained by 
the Doppler shifts of neutral hydrogen, is considered to 
be sampling a component of z-velocity if the mushroom- 
shaped structure is tilted relative to the line of sight. In or- 
der to show the velocity structure from our simulation, we 




Fig. 3. Averaged velocity across the cap as a function 
of R in the case of v zc0 = 10c s0 ~ 100 km/s. The upper 
panel shows the velocity along z — 5.0H and the bottom 
panel along z — 3.0H - 



calculate the averaged z-component of the velocity along 
a line of sight perpendicular to the R — z plane. The av- 
eraged velocity is calculated in the following manner: 



where I is the length along a line perpendicular to the 
R — z plane. The integration was also done only for the 
cold gas in the same way that we get the column density. 
Figure 3 shows the averaged velocity across the cap. The 
top panel shows the velocity along z = B.OHq, and the 
bottom along z = S.OHq. At the lower portion of the cap 
(z = 3.0-ffo), the flow goes up near the axis but goes down 
near the edge of the cap, while all of the cap is gradually 
going up in the upper portion (z = H.OHq). The velocity of 
the upper portion (v z ~ 0.4c s ) is in between the maximum 
and minimum velocities of the lower portion of the cap. 
These velocity structures are consistent with Figures lc 
and Id of English et al. (2000). 

3.2. The case ofv ZC Q = 5c s q — 50 km/s 

Here, we will show the case of v zc q — 5c s o ~ 50 km/s. 
The time evolution of the cloud is similar to the case of 
v Z co — 10c s0 ~ 100 km/s. Figure 4 shows the snapshot 
of density and column density at t = 2.75io- We can see 
the similar mushroom-shaped structure around t — 2.75to, 
though the size of the mushroom is a little smaller than the 
case of v zc o = 10c s o — 100 km/s. The maximum height 
of the mushroom at t = 2.75t is about 3H — 420 pc, 
which is almost the same size as that of GW 123.4-15. 
The width ratio of the cap and stem is about 2:1 and the 
mass ratio is about 4:1, which also show good agreements 
of the observation. 

Figure 5 shows that averaged velocity across the cap 
at t = 2.75t . The top panel shows the velocity along 
z = 2.0Ho, and the bottom along z = lAHo- The lower 
portion of the cap shows that the flow goes up near the 
axis but goes down near the edge of the cap, while the 
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Fig. 4. Results in the case of v zc q = 5c s q ~ 50 km/s. Left: density at t = 2.7bt$. Right: column density of low 
temperature gas (T < 0.87i) at t — 2.75t . The density is normalized by p the initial density at z = 0. The column 
density is normalized by PqHq. 



upper portion of the cap shows the flow goes up in wider 
region than that of the lower cap. This is also consistent 
with the observation. Movies of the density and column 
density for this model further illustrate the evolution. 

3.3. Parameter study 

We have performed a parameter study for the velocity 
(vzco) and density (p c0 ) of the impact cloud, as well as 
the specific heat ratio (7). 

When the velocity is 3c s0 (— 30 km/s), the mushroom- 
shaped structure does not appear. The cloud docs not 
reach enough height after it goes though the dense sheet, 
so that the cap and stem structure is not apparent (the 
cap is prominent, but the stem is not). Therefore, we need 
at least about 50 km/s to get a mushroom-shaped struc- 
ture when the density of the cloud is the same order of 
that of the dense sheet on the Galactic plane. 

If the density of the cloud is much less than that of the 
Galactic plane (p c < O.lpo), our simulations show that the 
cloud does not go through the plane even in the case of 
v Z cO — 10c s o — 100 km/s. Therefore, the cloud density 
should be of the same order as that of the dense sheet 
of the Galactic plane unless the velocity of the cloud is 
extremely high. 

In order to see the effect of cooling mechanisms (ra- 
diative cooling and/or dissipation of turbulence), we ran 
simulations with 7 = 1.05. Even for this small specific 
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Fig. 5. Averaged velocity across the cap as a function of 
R in the case of v zc q — 5c s q ~ 50 km/s. The upper panel 
shows the velocity along z = 2.0H n and the bottom panel 
along z = IAHq. 



heat ratio, we had similar mushroom-shaped structures 
both for high and low velocity cases. Therefore, we think 
that the effect of cooling does not affect the basic result. 
However, we found some differences: (1) the impact cloud 
reaches a higher position (about 20% higher) than that of 
the adiabatic case after it goes through the dense sheet. 
(2) Both stem and cap become a little narrower (by about 
10%) than those of the adiabatic case. 
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4. Discussion and Conclusion 

Based on our simulations, we conclude that a cloud colli- 
sion with the Galactic disk is one of the possible models of 
the mushroom-shaped structure, GW 123.4-1.5. The den- 
sity and column density clearly show a cap and stem struc- 
ture which resembles the observed structure for both cases 
of a high- velocity cloud (v zc o = 10c s o — 100 km/s) and an 
intermediate-velocity cloud {v zc q — 5c s q — 50 km/s). The 
case of v zc0 — 5c s0 ~ 50 km/s shows a better agreement 
with the size of the mushroom-shaped structure as well as 
the mass ratio of the cap and stem. It may mean that GW 
123.4-1.5 is created by the collision of an intermediate- 
velocity cloud into the Galactic plane, though the stem in 
our simulation is a little thicker than that of the observa- 
tion. 

However, mushroom-shaped structures may not be 
commonly observed unless the densities of the impact 
cloud are comparable to the density of the Galactic 
plane and the frequency of the impacts is high. In or- 
der to get a mushroom-shaped structure, we found that 
the cloud density should be at least the same order 
as the disk density. Given a disk number density ~ 1 
cm~ 3 , the column density of the impact cloud is > 10 20 
cm -2 . Observations show that such large column densities 
are not common both for high-velocity clouds (Wakkcr 
& von Woerden 1997) and intermediate- velocity clouds 
(Benjamin & Danly 1997). Nevertheless, there are some 
exceptions, such as the Ursa Major intermediate-velocity 
cloud (Snowden et al. 1994; Benjamin et al. 1996) which 
has velocity ~ 50 km/s, column density ~ 2.0 x 10 20 cm~ 2 , 
and estimated size (~ 15 x 50 pc) which implies a number 
density ~ 1 cm~ 3 . The numerical simulations also show 
that the lifetime of the mushroom-shaped structure is only 
about a dynamical time, ~ 10 7 years. These results may 
be some of the reasons that only one mushroom-shaped 
structure has yet been found near the Galactic plane. 

The cloud collision model in this paper is an alternative 
to the buoyant model. The velocity structures across the 
cap in this model are consistent with observations of GW 
123.4-1.5, but it is not obvious that the buoyant model 
has a similar velocity structure. More detailed kinemat- 
ical analysis of GW 123.4-1.5 would distinguish the two 
models. In addition to the velocity, the stem near the mid- 
plane of the Galaxy is clear in our model, while it is not 
clear in the buoyant model. Our model shows that the 
stem directly connects to the dense gas near the midplane 
of the Galaxy, which seems to be the same as the obser- 
vation. Moreover, our model predicts small structures on 
the opposite side of the Galactic plane (in Figure 2 around 
z = — 1.5.ffo) an d in Figure 4 around z = ~1.0H ). These 
structures may not be significant in current observations, 
but may be an interesting target in the future in order to 
distinguish the two models. 

We have assumed a Galactic atmosphere with two ef- 
fective temperatures. The resulting dense sheet of gas near 
the midplane contributes to the formation of the stem of 
the mushroom-shaped structure. In our model, the stem 



originates from material lifted up from the dense sheet, 
while the cap is formed primarily by the breakup of the 
impact cloud. A dense stem is actually not clear in single- 
temperature models that we performed with the same 
initial impact clouds. A distinct transition between cold 
(weakly turbulent) and warm (highly turbulent) gases is 
not obvious from observations which measure the inte- 
grated gas along the line of sight. However, the mushroom- 
shaped structure we have obtained in our simulations may 
suggest that the cold gas clouds which are located near the 
midplane of the Galaxy are locally well separated from the 
warm gas. 

Future work can improve upon this model by including 
a self-consistent treatment of turbulent pressure, and the 
effect of radiative cooling of the thermal gas. Moreover, 
the magnetic pressure is comparable to the effective tur- 
bulent pressure in the ISM, and a collision exactly perpen- 
dicular to the Galactic plane would be rare in a real situ- 
ation. Three-dimensional simulations including the mag- 
netic field, heating and cooling, and turbulent gas would 
be an ultimate work for the future in order to have more 
detailed comparison between the numerical simulations 
and observations. 
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